Spherically symmetric massive scalar fields in GR 
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First we review some of the attempts made to find exact spherically symmetric solutions of 
Einstein field equations in the presence of scalar fields .Wyman solution in both static and non 
static scalar field is discussed briefly and it is show that why in the case of non static homogenous 
matter field, static metric can not be represented in terms of elementary functions. We mention 
here that if our spacetime be static, according to EFE there is two option for choose scalar field 
■ matter: static (time independent) and non static (time dependent). All these solutions are limited 

D ' to the minimally coupled massless scalar fields and also in the absence of the cosmological constant 

. Then we show that if we are interesting to have a homogenous isotropic scalar field matter one can 
^ — . construct a series solution in terms of scalar fields mass and cosmological constant. This metric is 

static and posses a locally flat case as a special chooses of mass of scalar field and can be interpreted 
i i as an effective vacuum. Therefore mass of scalar field eliminates any locally gravitational effect as 

tidal forces. Finally we describe why this system is unstable in the language of dynamical systems. 
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I I. BACKGROUND 

> : 

Static spherically symmetric spacetimes are the simplest ones to study. Bcrgmann and Leipnik [1] are the first who 
studied solvability of field equations in the presence of scalar fields. They were not able to find an explicit solution due 
to their choice of a Schwarzschild like coordinate. Buchdahl in a series of papers [24] introduced and developed the 
reciprocal metrics idea as a generating method to construct new vacuum or non vacuum static solutions as a simple 
transformation of a known vacuum or non vacuum one. The higher dimensional extension of Buchdahl family has 
been introduced by Xanthopolous and Zannias [5]. Later Wyman [6] attacked this problem in a general framework 
and showed that apart from the Buchdahl family of solutions there is another family of static spherically symmetric 
solutions in which the scalar field is not static. Later critical phenomena in gravitational collapse were discovered by 
Choptuik [1315] in the model of a spherically symmetric, massless scalar field minimally coupled to general relativity 
. The scalar field matter is both simple, and acts as a toy model in spherical symmetry for the effects of gravitational 
radiation. Important analytical studies of gravitational collapse in this model have been carried out by Christodoulou 
[1622]. For other numerical work on this model, see [2327]. In this work after a brief review of some of the previous 
studies we show that there is no closed form of Wyman class and then in section III first we derive EFE for a massive 
minimally coupled scalar field in the presence of cosmological constant in a spherically symmetric space time. Then 
by supposing that our scalar field has small time variations such that we can disregard its time dependency, and 
corresponding to it a static metric, we derive a series solution for our model. We mention here that if we assume that 
our field remains static, pressure and energy density describe a stiff matter and our spacetime posses homogeneity . 
This series solution has positive spatial parity and by choosing a suitable mass of scalar field is locally flat. Finally 
we study the stability of these solutions as a dynamical system and as an important result we show that Wyman 
minimally coupled solutions is unstable in the language of dynamical systems. 

II.NON INTEGRABILITY OF NON STATIC SCALAR FIELDS IN SPHERICALLY SYMMETRIC 

SPACETIMES 

We start from the Lagrangian of a static massless minimally coupled spherically symmetric scalar field with action: 

S= f (-g)HR + ^^Md 4 x (1) 
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where in it (/i, v = 0, 1, 2, 3). We can write the static spherically symmetric solution in schwarzschild gauge as: 

ds 2 = -e u dt 2 + e x dr 2 + r 2 {dd 2 + sin 2 $d<p 2 ) (2) 
Einstein and scalar field equations are: 

□0 = (3) 

Buchdahl in series of papers introduced a new transformation which we called reciprocal transformation. Im- 

posing Buchdahl transformations to the Schwarzschild as the seed metric ,we end up with a one parameter family 
of exact solution for a minimally coupled massless static scalar field . As another example of the same technique a 
static solution with axial symmetry in weyl coordinates can be obtained [5[ starting from the Weyl solution. The one 
parameter Buchdahl solution in spherically symmetric coordinate is given by: 

ds * = _ (1 _ ™)Ptf + (1 - —yVdr 2 + r 2 (d$ 2 + sinHdAp 2 ) 
r r 
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Field equations are now simplified to; 

Run = -^-^w (6) 

Wyman showed that Q in Buchdahl formalism one family of solutions is missing and Yilmaz [3], Szekeres Q and 
Buchdahl solutions are all sub classes of his static family. Here we restrict ourselves to Wyman's non static spherically 
symmetric scalar field family of solution ; i.e the case in which A = A(r), v = v(r) , 4> = 4>{t) and the metric functions 
(2) satisfy the following differential equation equations: 

v + A = ure x ~ v 
z/-A = 2(l-e A )/r 

Eliminating v we find the following for A; 

A +^--(e A -l) + ^(e A -l)( e A -2) = (8) 

Wyman has stated that the above equation (8) can not be integrated and by imposing suitable boundary conditions 
derived a series solution for the metric functions. Now using results from differential equations we show that this 
family could not be written in closed form. Using the following change of variables : 

t ■ dX \ 
r = e , u — A = — , x — e , (9) 

Equation (8) is transformed into: 

du 4 — 3a; 2 , . . 

u— — = u 0-l)(x-2) (10) 

dx xx 

This differential equation is of 2-nd Abel type (class A) [9] . It is well known that for a general equation of the form: 

(y(x) + g(x))^- = f 2 (x)y(x) 2 + h(x)y(x) + f (x) (11) 
dx 

There are solutions only in the following two special cases: 
Case(l): If 

/o(a0 = fi{x)g(x) - f 2 (x)g(x) 2 

y(x) = -g(x) , 12 x 

y(x) = e/A(*)«te( Cx + f (-f 2 ( x )g(x) + h(x))e- f f*^dx) 
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Case(2): If 

/iO) = 2f 2 (x)g(x) -g{x)g{x) 
a{x) =e -zfMx)dx 

P(xf = a { xf 9{ x? + 2a(x) Cl + 2a(x) J ^ £f£ )dx)2 («) 

a(x) 5 (x) + A//3(a;) 

= , x 

Since Wyman equation does not belong to any of these cases , one can not derive an exact closed form for it. 

In this section we obtain a new series solution for a massive-minimally coupled scalar field in the presence of a 
cosmological constant. According to Wyman we take our spherically symmetric static spacetime metric in Schwarzschild 
gauge, and require the metric functions to be regular near the origin . So that the scalar field and metric functions 
have a Taylor expansion near the origin. The field equation are now given by: 

Rfj,u ~ Ti^-Siiv ~ ^-9nv = ~kT^ u 

2 1 (14) 

KT\u> = ~lA<l>;n4>;i> - -^9iiv{4>, a 4> a ~ m 2 (j) 2 )} 

Where we have introduced a parameter fi in the energy-momentum tensor to make sure one that our solution reduces 
to the vacuum solution in the absence of the scalar field, we take it 1 for simplicity. The equivalent equation can be 
written as: 

R» v = -ASV _ 4>-^4>-m + ^m 2 <j) 2 g^ (15) 
Assuming that our spherically symmetric metric in the the Schwarzschild gauge is of the general form: 

ds 2 = e u dt 2 - e x dr 2 - r 2 (d6 2 + sin 2 6d<p 2 ) (16) 
The scalar field equation of motion is: 

nd} + m 2 (j) = 0, (17) 
and equation (15) reduces to the following differential equations: 

-^-e~ x (2rv" + r(y) 2 - rv X + Av) = -A + \-m 2 d) 2 



— {2rv + r{v') 2 - rvX' - 4A') = Ae A - 2 - -m 2 (f> 2 e x 

^e~ x (-rX - 2e A + rv + 2) = -Ar 2 - ^m 2 <j) 2 r 2 
d 2 6 „ .1 v — X .dd 9X , 



(18) 



From (18) after some little algebra we find: 



v + X = -r6 2 



■ 2 1 (19) 

v - X = -[e x {\ - Ar 2 ) - 1 - ^m 2 d> 2 r 2 e x ] 

There is no simple reduction of this system to a one variable solvable differential equation. Actually field equation 
can be writhen in the following complicated non linear form: 

^ 2 _ h 3 (f> + 3h 2 <p (20) 
ho 
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in which: 

ho = rcf>" (rm 2 cf) + <f> (-2 + 2Ar 2 + m 2 r 2 (j) 2 )) 

h 3 = (-2 + 2Ar 2 + m 2 r 2 <j) 2 )(j)" + ^-m 2 r*(j)(f> + (m 2 (j) 2 + ^— + 2A) + \m 2 4> (21) 

h 2 = (4 - 4Ar 2 - 2m 2 r 2 ^ 2 )^' - 2m 2 r^ 

It does not seem one could be able solve (20) analytically, so we try a semi-analytic solution by imposing the same 
simple boundary conditions on metric as Wyman did (f|: 

1/(0) = 0,A(0) = 0,0(0) = qj(0) = (22) 

These boundary conditions are chosen such that we have a regular solution at the origin, q is interpreted as a new 
parameter of our two parameter (q, m) family of exact solutions, apart from the mass of the scalar field and the 
cosmological constant. Applying the above boundary conditions in equation (19) and taking successive derivatives 
of equations and evaluating them at the origin we can write the following series for the scalar field and the metric 
functions 12j : 

<P(r) = q(l+ Yl M?V) 

n=2,4,6,... 

e" = l- £ b n (q 2 )r n (23) 

n=2,4,6,... 

e A = l+ J2 C «(</V 

n=2,4,6,... 

We note that only even values of n contribute in the solution. This means that our functions remain unchanged 
under parity transformation r — > — r. We can conclude that our metric functions can be written in general form 
X(r) = 1 ± f(m 2 q 2 ,Ar 2 ). The first few values of the coefficients in (23) are give below: 

a 2 (x) = ^m 2 

a 4 (aO = ^ (^m 4 x + ^m 4 + ^Am 2 ) 

b 2 (x) = hm 2 x + 2A] 

6 ! ; (24) 

c 2 (x) = hm 2 x + 2A] 
6 

. . l r 16 2 2 42 8 2 a 8 n, 
c±(x) = — 1^5 x + 3™ x 3 TO 3 ' 

If q = the solution reduces to that of the series expansion for de-Sitter solution . It is clear that there is no consistent 
Wyman solution of 2-nd class. Calculating all the Rieman Tensor components we observe that the only non zero 
component up to 4-th order in r (near the origin) is given by: 

R% = ^m 2 q 2 - iA (25) 



So if the mass of the scalar field satisfies m = q 1 V2A this component vanishes and the spacetime becomes fiat and 
every local effect of gravity as tidal forces vanishes. 



MASSIVE SCALAR FIELD AS A DYNAMICAL SYSTEM 



In this section we apply a powerful method to analyze the stability of the scalar field solutions obtained in the 
previous section . This method for an autonomous ODE has already been employed by Coley f5~oj ] in the context of 
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cosmology. To do so we introduce some concepts from nonlinear complex systems ll| . In the language of dynamical 
systems our differential equation could be represented as follows : 



x = f{t,x) 
f : [0, oo) x D — > K" 
D = {xe M" HI x || a < 0} 

We also assume that x = be the equilibrium point of the system at t = that is: 

/(t,Q) = 0,Vt>0 

• rail 



(26) 



(27) 



We assume that the Jacobian matrix ; [-^] is bounded with respect to t on the set D and smoothly satisfies Lipschit'z 
lemma , that is: 



f(t,x) - f(t,y) \\< L || x- y ||, || x \\ p = (]T | x t \ p )p,1< P < oo 



(28) 



One can state the following theorems about this system. 

Theorem (1): In the neighborhood of any stability point the system can be approximated by a linear one i.e;: 

x = A(t)x (29) 
Theorem (2): Assume that x = be an equilibrium point of the nonlinear system x = f(t,x) that is : 

/:[0,oo)xD — > K" 



D = {xeR n HI x \\ 2 <r}, 



(30) 



and define the matrix; A(t) 



df(t,x) 



\x=o ■ If the origin is an exponential stable equilibrium point of the linear system 



(29) then this point is the exponential stable equilibrium point of the non-linear system (26). 

Now we apply these theorems to our system. First we rewrite held equations through the following variables:, 



Xi = v(r),x 2 = \{r),x 3 = 4>(r),x 4 



d<t>{r) 
dr 



(31) 



as follows 



e X2 - 1 



Xi 



+ Are X2 r{x\ + m 2 e X2 xl) 



e X2 — 1 1 
x 2 = Are X2 r{x\ — m 2 e X2 xl) 



(32) 



X4 



-e X2 + 1 1 

Are X2 + -m 2 rx 2 e X2 ]xi - m 2 x 3 e 



X 3 = X4 



Equilibrium point of the above system lies at: 

x\ = Xi,X2 = — ln(l + Ar 2 ), x 3 = 0, X4 = 



(33) 



This point of 4-dimensional phase space corresponds to a t = cte hypersurface with constant negative intrinsic 
curvature k = —A with metric: 



ds 



dr 2 



t=cte l + A? , 2 



+ r 2 (d9 2 + sin 2 6d(p 2 ) 



(34) 



We assume that the equilibrium point (solution) is perturbed as x.- L = x\ + Uj. The Jacobian matrix can be easily 
calculated and the result is: 



At 



dx-i 





' r- 1 





" 




r- 1 








(t,x)] = 














1 




. 


m 2 
1+Ar 2 


_ 



(35) 
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We can find an exact solution for perturbations of our dynamical system variables (29) according to the following 
system of ordinary differential equations, using (29) ; 



1 

Ul = -u 2 
r 

1 

U 2 = -Ul 

r 

u 3 = u 4 



1 + Ar 2 



(36) 



From these set of equations we immediately obtain: 



Mi = c x r + c 2 
u 2 = c\r 

u 3 = (l + Ar 2 )(c 1 F([a 1 ,a 2 ],l/2,-Ar 2 )+ 
c 2 F([ ai + 1/2, a 2 + 1/2], 1/2, -Ar 2 )) 

d,U4 

u A = — — 

ar 



(37) 



In which: 



° z k — — 



F(n, d,z) = J2jj Y[Mk ]]_(dj)k 

k=0 ' i=l j=l 

T(k + d 3 ) 

n = [ni,n 2 , ...],p = nops(n),d = [di,d\, ...],q = nops(d) (38) 

3\/A + iV4m 2 - A 



ai = 1/4 
a 2 = -1/4 



^/A 

3\/A + iV4m 2 - A 



and the symbol nops means the number of operands of an expression. If some is a non-positive integer, the series 
is finite (that is, F(n, d, z) is a polynomial in z). If some dj is a non-positive integer, the function is undefined for all 
non-zero z, unless there is also a negative upper parameter of smaller absolute value, in which case the previous rule 
applies. For the remainder of this description, assume no rii or dj is a non-positive integer. When p < q , this series 
converges for all complex z, and hence defines F(n, d, z) everywhere. When p = q+ 1 , the series converges for \z\ < 1 
. F(n,d,z) is then defined for \z\ > 1 by analytic continuation. The point z = 1 is a branch point, and the interval 
(1, oo) is the branch cut. When p > q + 1 the series diverges for all z ^ . In this case, the series is interpreted as 
the asymptotic expansion of F(n, d, z) around z = . The positive real axis is the branch cut. 
Now the main result comes from the fact that the solution (37) is not asymptotically convergent i.c; 

Ui(r -4- oo) = oo (39) 

So although we treated the system in the linear approximation, we deduce that it is not stable . Therefor it is 
guaranteed that our nonlinear system is not stable too. For a massless (Wyman) solution in the absence of the 
cosmological constant we can refer to the system of equations (36) and its solutions (37). Obviously in this case the 
perturbations diverge and so this family is not stable as well. 
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